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The so-called extended law of corresponding states, as proposed by Noro and Frenkel 
[J. Chem. Phys. 113, 2941 (2000)], involves a mapping of the phase behaviors of sys¬ 
tems with short-range attractive interactions. While it has already extensively been 
applied to various model potentials, here we test its applicability to protein solutions 
with their complex interactions. We successfully map their experimentally deter¬ 
mined metastable gas-liquid binodals, as available in the literature, to the binodals 
of short-range square-well fluids, as determined by previous as well as new Monte 
Carlo simulations. This is achieved by representing the binodals as a function of the 
temperature scaled with the critical temperature (or as a function of the reduced 
second virial coefficient) and the concentration scaled by the cube of an effective 
particle diameter, where the scalings take into account the attractive and repulsive 
contributions to the interaction potential, respectively. The scaled binodals of the 
protein solutions coincide with simulation data of the adhesive hard-sphere fluid. 
Furthermore, once the repulsive contributions are taken into account by the effective 
particle diameter, the temperature dependence of the reduced second virial coeffi¬ 
cients follows a master curve that corresponds to a linear temperature dependence 
of the depth of the square-well potential. We moreover demonstrate that, based on 
this approach and cloud-point measurements only, second virial coefficients can be 
estimated, which we show to agree with values determined by light scattering or by 
DLVO-based calculations. 
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I. INTRODUCTION 


Proteins show a rich phase behavior; they can crystallize,^ undergo liquid-liquid phase 
separation (LLPS),®' or form aggregate^®' or hbrils.®^ It depends on the protein-protein 
interactions which can be tuned through the solution conditions, namely thermodynamic 
parameters, such as temperature, pH, salt and protein concentration.'^®!' Controlling the 
phase behavior of protein solutions is of fundamental importance in medicine, materials 
science, physics and several other disciplines, as illustrated by the following examples: Con¬ 
densation of protein phases, as observed in LLPS and hbril formation, plays a major role in 
the pathogenesis of several human diseases,'!®!' such as cataract formation and Alzheimer’s 
disease.'!®^ Unwanted protein phase transitions can severely impede drug delivery and purih- 
cation in pharmaceutical treatments and food engineering.'!®®^ Moreover, conditions under 
which high quality protein crystals grow, which are an essential prerequisite for diffraction 
studies, are related to the protein phase behavior.'®' 

Protein solutions represent a highly complex system comprising protein molecules, water 
and usually also salts and further additives. Hence, there are a number of different contribu¬ 
tions to intermolecular interactions among proteins, such as electrostatic interactions, van 
der Waals forces, hydrophobic interactions and hydration, ion-dispersion forces, and ionic or 
other bridges, etc., which are often highly directional. Due to this complexity, the protein- 
protein interaction potential U{f), where r*is the center-to-center vector of protein molecules, 
is usually not known in detail. Nevertheless, various models from liquid-state theory and col¬ 
loid science, in which protein solutions are represented as effective one-component systems, 
have successfully been applied to study protein phase behavior. In many cases, protein- 
protein interactions are modelled by a hard-core repulsion, a short-range attraction and a 
long-range repulsion. For example, the Derjaguin-Landau-Verwey-Overbeek (DLVO) poten¬ 
tial consists of a hard core, a short-range attractive van der Waals contribution and screened 
Coulomb interactions, thus directly taking into account the charges of proteins and charge 
screening due to salts. Several authors have studied protein phase behavior using DLVO- 
based models.'®® Applying modihed DLVO theory, salt concentration and pH effect^®® 
as well as the effect of glyceroP®!' on protein crystallization could be explained. Although 
DLVO is not able to explain, e.g., salt-specihc effects,'®'®' it has been arguecp!' that the 
functional form of the (attractive part of) DLVO potential is flexible enough to also account 
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for non-DLVO forces, such as hydration, if the interaction parameters (e.g., the Hamaker 
constant) are suitably adjusted. Especially in Monte Carlo simulations, the hard-core at¬ 
tractive Yukawa and the Asakura-Oosawa (AO) potentiaP^^ have also been 

applied to model protein-protein interactions in the presence of additives. The simplest 
model accounting for attractive forces is the square-well (SW) fluid, short-range SW poten¬ 
tials have frequently been used to model protein-protein interactionsincluding the 
analysis of small-angle scattering data of protein solutions.!^ Moreover, Baxter’s model, i.e., 
the adhesive hard-sphere (AHS) potential, which represents a specihc limit of the SW poten¬ 
tial, has also been applied to study protein-protein interactions.®^^® Recently, simulation 
work on anisotropic interactions of patchy particles accounting for the high directionality 
of protein interactions was successfully compared with binodals of protein solutions.®®® 
Although the effective interaction modeld®® involve a (crude) simplihcation of the interac¬ 
tions in protein solutions, reasonable modelling of protein phase behavior is achieved if the 
short-range nature of the interactions is considered. 

An integral parameter, which is independent of the specihc shape of the interaction po¬ 
tential, is the second virial coefficient B 2 . It provides a quantitative measure of interactions: 
positive (negative) values of B 2 indicate net repulsive (attractive) interactions. For interac¬ 
tion potentials with spherical symmetry, 17(r), the second virial coefficient is given by 

U{r 


B2 = 271 


1 — exp 


^rT 


r dr. 


( 1 ) 


where /cr is Boltzmann’s constant and T the absolute temperature. The normalized second 
virial coefficient 62 is dehned as 

B2 


bo = 




( 2 ) 


where B^^ = |vrcr^ is the second virial coefficient of a hard sphere system with particle 
diameter a. Second virial coefficients can be determined experimentally, e.g., by scattering 
techniques and osmometry.®®®® 

Some empirical “laws” linking the phase behavior of colloids and proteins to the second 
virial coefficient have been established. Based on experiments, George and Wilson® have 
shown that optimum solution conditions for the crystallization of a number of proteins are 
characterized by a range of 62 values, the so-called crystallization slot®® 


- 10 < bo < -1 


r\j ^ r\j 


(3) 
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For lower 62 values, attractions between protein molecules are so strong, that amorphous pro¬ 
tein aggregation occurs, while, for larger 62 values, repulsions dominate and crystallization 
is too slow or impossible. 

For fluids with short-range attractions, quantitative scaling laws for the gas-liquid binodal 
near the critical point have been proposed based on simulation results. Vliegenthart and 
Lekkerkerker (VL)P51have found that, for many model potentials, the second virial coefficient 
takes an approximately constant value at the critical temperature T^. 


h(T,) = 6 ^ « - 1 . 5 . 


( 4 ) 


Moreover, Eq. Q can be used as a predictor for the critical temperature: if, for a given 
temperature T, one hnds h 2 {T) —1.5, then Eq. (|^ implies T ^ Tc. Recently, Wolf et 

al.^ conhrmed the VL Ending. They argue that the osmotic pressure of a protein solution 
close to its critical point is low and hence it can be expanded into a virial series truncated 
after third order in density, yielding an approximate relationship between the critical volume 
fraction 0 c and second virial coefficient at the critical temperature: 


b 


C 

2 



(5) 


which, based on Eq. (|^, suggests 0c ~ 0.17 in agreement with other findings (e.g.. Sec. 

Fig. 1^;). 

Applying the Weeks-Chandler-Andersen (WCA) method,l^^^ according to which an in¬ 
teraction potential U{r) = f/rep(’") + f^attr(’") is separated into a contribution containing all 
repulsive interactions, Urepir), and one containing all the attractive parts of the potential, 
Ua.ttri'f'), Noro and Frenkel (NF)P^ have suggested that many model potentials with suffi¬ 
ciently short-range attractions can be characterized by only three quantities: (i) an effective 
hard-core diameter Ueff taking into account all repulsive interactions between the particles,^ 
which is given by the relationship 


Ceff — 


1 — exp 


Urep (^) 

knT 


dr. 


( 6 ) 


(ii) an energy scale e being a measure for the strength of the attraction, and (iii) the reduced 
second virial coefficient b^, which reads 


6 * = 


B2 


Bp 






(7) 
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It is related to the range (and strength) of the attractions only, becanse the contribntions of 
the repnlsive interactions to B 2 are balanced by considering cXeff instead of a. Based on their 
mapping onto an eqnivalent sqnare-well flnid, Noro and Frenkel proposed an extended law 
of corresponding states (ELCS): Many colloidal systems with the same rednced temperatnre 
T* = k^T/e, density p* = pa^f^ (where p is the particle nnmber density), and second virial 
coefficient obey the same eqnation of state, snch that the compressibility factor z is a 
fnnction of only three parameters: 


n 


2; = 


pkjiT 


= z{T\p\h%), 


( 8 ) 


where 11 is the osmotic pressnre of the system. 

While for many model potentials, inclnding HCAY, AO and patchy particles, the 
applicability, possible extensions and restrictions of the ELCS have been stndied by 
simnlations,^21IiniSSEI]^Pg ELCS has hardly been assessed for systems as complex as proteins.^ 
Thongh it has been fonnd that scaling gas-liqnid binodals of protein solntions to the critical 
points resnlts in a master cnrve,®E2lniost experimental valnes of the second virial coefficient 
62 are systematically smaller than those predicted by VL and knowl¬ 

edge, systematic attempts to clarify the differences between the qnantitative predictions 
of the ELCS and experimental data of the gas-liqnid binodal of protein solntions and the 
corresponding temperatnre dependent second virial coefficients have not been nndertaken 
yet. 

In this work, we aim for a direct and qnantitative comparison of the predictions of the 
ELCS with the experimentally observed phase behavior of protein solntions, namely the 
metastable gas-liqnid binodal, and the temperatnre dependent second virial coefficient. This 
comparison is based on experimental data available in the literatnre that provide both the 
phase behavior and second virial coefficients and also on previons and new Monte Carlo (MC) 
data on short-range SW flnids. We start (Sec. 0 by applying the ELCS to short-range SW, 
AHS, HCAY and AO flnids. In the &2 plane (where (j) is the particle volnme fraction), the 


binodals of all these flnids match. Then, in Sec. Ill, the experimentally determined binodals 
of protein solntions are compared with those of SW flnids. Inclnding the repnlsive part of the 
protein-protein interactions throngh an effective hard-core diameter aes yields a reasonable 
collapse of the experimental protein data onto a single binodal in the T /Tc — cpes plane (where 


(j)es is the effective volnme fraction based on (Jeff). In Sec. IV, we snbseqnently investigate the 
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temperature dependence of the second virial coefficient. If the effective hard-core diameter 
(Teff is taken into account, for a variety of solution conditions close to the critical temperature 
Tc, the experimental values of the reduced second virial coefficient b^iT/T^) are found to lie 
near a master curve, that corresponds to a linear temperature dependence of the depth of the 
SW potential. In Sec.|V| then the b^iT/Tc) values are used to plot the experimental protein 
binodals in the b^ — (pes plane. They match with the binodals of the short-range SW and the 
AHS fluids, thus conhrming the applicability of the ELCS to protein solutions. Furthermore, 


in Sec. VI we show that, based on the DLVO potential, the 62 values of protein solutions 
can successfully be calculated for different salt concentrations and pH values. However, 
using DLVO theory the effective hard-core diameter is underestimated, indicating that the 
protein-protein potential contains repulsive contributions beyond the electrostatic repulsion. 


We hnally demonstrate in Sec. |VH| how our approach can be exploited to determine 62 values 
based on cloud-point measurements only, which are shown to quantitatively agree with values 
obtained by light scattering experiments or by calculations based on the DLVO potential. 


We conclude by summarizing our Endings in Sec. VHI 


II. SQUARE-WELL FLUIDS AND THE ELCS FOR SIMPLE FLUIDS 

The square-well (SW) potential has a simple mathematical structure and provides an 
unambiguous definition of the interaction range. It is often used as a simple model for the 
interactions governing the phase behavior of protein and colloid systems. In order to provide 
a basis for the mapping of the protein phase behavior onto an equivalent SW potential, 
following Noro and Frenkel,^ we discuss features of short-range SW fluids in this section. 
This discussion is based on our new MC simulation data as well as results available in the 
literature. 

The coexistence of two fluid phases, a gas phase and a liquid phase, in a one-component 
system is only possible if attractive particle-particle interactions are present. The simplest 
model showing this scenario is probably given by the square-well (SW) potential Usw{r). 
It consists of a hard-core repulsion of range a (the diameter of the particle), which leads 
to excluded volume effects, and a constant attractive part, which extends to a distance Xa 
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from the center. Thus, the SW potential is mathematically dehned as 


oo, 

Usw{r) = -e, 

0 , 


r < a, 
a <r < Act, 
r > Act, 


( 9 ) 


where e > 0 is the depth of the square well, quantifying the strength of attraction.® (Note 
that the dimensionless width of the SW, 5 = A — 1, is also used in the literature.) 

The phase behavior of square-well fluids is governed by the ratio of the range Act and 
the particle size ex.® It has been shown®®® that for A < 1.25 gas-liquid phase separation 


is metastable with respect to the fluid-solid equilibrium. We will refer to SW fluids with 
A < 1.25 as short-range SW fluids. Simulations indicate that also for a Lennard-Jones 
potential metastable gas-liquid phase separation occurs for similar values of the range.® 
Experiments and theory on colloid-polymer mixture^®® have shown that the topology of 
the phase diagram changes when the range of interaction is about A ~ 1.25. Moreover, 
the limit of short ranges A is relevant for exploring the relation of the SW potential to the 
adhesive hard-sphere (AHS) potential,®®® which is defined by an infinitesimal interaction 
range A in combination with an infinite interaction strength such that 62 remains finite. 

In Figure [T| we present new MC simulation data for A = 1.05 and 1.1 as well as simulation 
and perturbation theory data of square-well fluids collected from the literature.®®®®®® 
Figure [pi shows gas-liquid coexistence curves of short-range SW fluids as obtained by MC 
simulations, where the critical points have been estimated using scaling laws and the law of 
rectilinear diameters.® It is important to note that data in the very vicinity of the critical 
point (0.95 < T/Tc < 1) is quite scarce. Duda!®was able to perform NVT—MC simulations 
of SW fluids with A = 1.05 and 1.10 in this temperature window. Duda’s values for the 
critical parameters are slightly off those obtained by Largo et ah® Thus, we have done 
simulations using the Gibbs ensemble technique,® as described previously.® MC runs for 
SW fluids with A = 1.05 and 1.1 were performed (Tab. 0 in order to establish the SW 
binodal in the vicinity of and to provide estimates for the critical parameters, which 
agree with previous result^ (Fig- [^)d). The error bars of the new simulation data (red 
symbols with black frames) are smaller than the symbol size used in the plots, except the 
one for the critical points. With increasing range A, the reduced coexistence temperature 
T* = k^T /e increases for all volume fractions studied due to enhanced attractions in fluids 
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of larger range. This trend is also observed in the critical temperatures T* (Fig. [^, inset). 
Furthermore, a slight shift of the critical volume fraction 0c to smaller values with increasing 
range A is observed (Fig. HI)- The width of the coexistence region, limited by the gas (low 
0) and the liquid (high 0) branches, slightly depends on A. 

Figure 03 shows the gas-liquid binodals of short-range SW fluids with temperatures 
normalized to the best estimates of the critical temperature. The gas branches of the 
different binodals seem almost to coincide within statistical error bars (in the low 0 regime), 
although the curves are not scaled to the critical volume fraction 0c. However, the liquid 
branch of the coexistence curve moves toward smaller 0 with increasing range A in this 
representation. 

Figure shows the best estimates of values of the volume fraction 0c at the critical 
point (black symbols) as a function of interaction range A, as determined by simulation and 
perturbation theory. We also include values reported for the AHS fluids, as obtained from 
the compressibility and energy route as well as simulation.^^^®® For a given range A, a 
variety of values of 0c (A0c ~ 0.02 — 0.05) have been reported for the SW fluid. 

It was recently proposed that the density at the critical point should be constant in the 
limit A —>■ 1 if measured in units of the average distance between two bonded particles.SSESl 
The simulation data of Figure 0:: have been scaled in this way (not shown). In the limit 
A —)■ 1, we observe an approximately constant value 

0c = 0c (^) = 0-288, (10) 


which can be regarded as an approximation to the critical volume fraction of the AHS fluid. 
Our value for 0° agrees with that reported by Largo et al.^ (0.289), but it is slightly larger 
than that for the AHS model predicted by simulations (0.266).^ For A > 1.5, the values of 
0c are systematically larger than expected based on Eq. (10) (Fig. 03, line). 

For a given reduced coexistence temperature T*, the normalized &2 value of a SW fluid 
of range A can be calculated analytically: 


b2{T*) 


1 - 0 - 1 ) 




( 11 ) 


Based on values for the reduced critical temperature T* (Fig. 03 , inset), we have calculated 
the second virial coefficient of SW fluids at the critical temperature = b 2 {T*) as a function 
of the range (Fig. W, black symbols). The spread in b^ for a given range A (A 62 ~ 1) reflects 
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Figure 1. Square-well (SW) fluids with depth e and range A (filled symbols), (a) Gas-liquid 
binodals of short-range SW fluids, shown as reduced temperature T* = k-^T/e vs. volume fraction 
(p. (b) Gas-liquid binodals of short-range SW fluids, re-plotted from (a), scaled to the critical 
temperature Tc. (Inset) Reduced critical temperature T* as a function of SW fluid range A. (c) 
Critical volume fraction cpc of SW fluids as a function of their range A. (d) Normalized second virial 
coefficient of SW fluids at the critical point, b^, as a function of their range A. In (c) and (d), open 
colored symbols represent values obtained in experiments on protein solutions (Tab. 
displayed at an arbitrary range. New MG data (red symbols with black frames) and data taken 
frp ^32 | 36|38 | 4u|tif|ti4 MlMM7H] as indicated. 


different values reported for the critical temperatures T* due to statistical and extrapolation 
errors as well as finite-size effects, which in general result in an overestimate of 

The recently proposed approximative relationl^ between and 0c, Eq. ([^, has been 
tested for the SW fluid. Based on the interpolation of simulation data of 0c(A) (Fig. [^, 
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Figure 2. Binodals of short-range square-well (SW) fluids in the b 2 — (p plane for 1.01 < A < 1.25 
(filled symbols, interpolated by a grey-shaded area as a guide to the eye) jSMQlZlHlS] addition, 
short-range hard-core attractive Yukawa (HCAY), Asakura-Oosawa (AO) and adhesive hard sphere 
(AHS) fluids are shown (open symbols, see legend) As suggested by the ELCS,^the data for 
these potentials can be mapped onto short-range SW fluids, i.e., onto the grey shaded Noro-Frenkel 
(NF) area. 

dashed line) and Eq. (|^, the valnes have been calculated (Fig. dashed line). Despite 
the simplicity of the argument, the values obtained in this way are even in semi-quantitative 
agreement with simulation data of the SW fluid, although they are slightly too large. Similar 
agreement has been found for data of short-range HCAY and AO potentials (not shown). 
Furthermore, to £t the data of Largo et al.^ (Fig. nil. stars), Gazzillo and Pini^have used 
a second order polynomial (Fig. Ef. solid line) 

^2(-^) = (^ 2)0 + Cl(A — 1) -|- C2(A — 1)^, (12) 

where (& 2 )o) Ci and C 2 are htting parameters to describe the range dependence of the T* 
(Fig.[g), inset, line). 

For the short-range SW fluids (Fig. &) , we have calculated binodals in the b 2 — cp 
plane, as shown in Figure i.e., we have replaced the T* axis with the 62 axis using 
Eq. ([^. In addition to the SW potential, model potentials of liquid state theory, such as 
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the hard-core attractive Yukawa (HCAY) potentiaP^®® and the Asakura-Oosawa (AO) 
potential,have also been used extensively to model interactions in protein solutions 
and colloidal suspensions. Here we focus on the gas-liquid phase separation only, although, 
for short-range potentials, crystallization intervenes and thus the gas-liquid coexistence be¬ 
comes metastable.^^^^^^®® The HCAY potential can be written as follows 


Chcay(?^) = 


oo, 


-ea 


exp [—K,(r—cr) 


r < (J, 
r > a, 


(13) 


where k quantihes the range of the HCAY potential. Sometimes, is used as an 

effective range of the HCAY model.^ The main quantity characterizing the AO potential, 
which is often used in the context of colloid-polymer mixtures, is the size ratio rj between the 
colloid and the polymer (for details see, e.g.j ^^ * ^^ ). Valadez-Perez et al.^^Slhave argued that 
short-range HCAY and AO potentials with kct > 3 and rj > 1.25, respectively, correspond to 
effective ranges Aeff < 1.25, assuming a constant 63 value of —1.5 (Eq. (|^). For the effective 


range, we obtain, based on Eq. (12), similar values (slightly larger ones for 77 , not shown). 

According to the Noro-Frenkel ELCS, it should be possible to map short-range HCAY 
and AO potentials onto short-range SW potentials, i.e., they should fall into the grey-shaded 
area of Figure which is a guide to the eye for the region covered by binodals of short-range 
SW fluids. For the parameter space of k and 7] stated above, we have calculated binodals in 
the b 2 — (j) plane (Fig. based on MC simulation data.^In addition, we include simulation 
results of the AHS potential.® All these data, i.e. all gas and liquid branches of the binodals, 
lie in the grey-shaded region apart from one exception. (The liquid branch for the AO fluid 
with T] = 1.25 lies slightly below the ELCS area; this value of r] seems to be too long-ranged.) 
The data shown in Figure [^provide further evidence for the mapping onto short-range SW 
fluids proposed by Noro and Frenkel.® Moreover, the values fall - by construction - into 
a fairly narrow range around —1.5, which is in line with the VL suggestion (Eq. 


III. LIQUID-LIQUID PHASE SEPARATION: GAS-LIQUID BINODALS 
OF PROTEIN SOLUTIONS 

The generic phase diagram of an aqueous protein solution contains a solubility (liquidus) 
line, indicating - for a given protein concentration Cp - the maximum temperature Txtai up to 
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which protein crystals are stable or - for a given temperature T - the maximum “solubility”, 
i.e. the crystallization boundary. For T > T^tai, a homogeneous protein solution without 
crystals represents the equilibrium state.^ (Note that thermal denaturation of proteins 
might interfere.) For temperatures below the crystallization boundary, metastable liquid- 
liquid phase separation (LLPS) of the protein solution into a protein-rich and a protein- 
poor phase is observed,® which is indicated by a dramatic clouding of the solution. At high 
protein concentrations, proteins might also form arrested gel or glassy states.®^^^^l^^® From 
a colloid physics point of view, proteins typically are regarded as colloids with short-range 
attractions and long-range electrostatic repulsion. In this respect, the LLPS coexistence 
curve corresponds to the gas-liquid binodal denoting the coexistence of a (dilute) colloidal 
gas and a (dense) colloidal liquid phase.^^ 

LLPS coexistence curves have been studied for proteins such as lysozyme, 7 -crystallins 
and antibodies.®^^®® In particular, lysozyme has often served as a model protein in 
physicochemical studies assessing the influence of different solution conditions, e.g., pH value, 
salt content or polymer concentration.^^HsziosHSl Moreover, lysozyme is one of the most studied 
proteins with respect to modelling of protein-protein interactionsj ^^ * ^^ * ^^ * ^^ * ^^ * ^^ ^ Therefore, 
we will focus on the binodals of lysozyme solutions here, although our analysis can also be 
applied to other protein systems. 

We have collected a number of gas-liquid binodals of lysozyme solutions for various 
solutions conditions from the literature.®2*^2*^^^^2ISSl Table we list solution conditions 
and interaction parameters of the experimental binodals, which we further analyze here. (We 
will use the values listed in Table for the DLVO model in Sec. |V^) The data sets cover a 
huge variety of different solutions conditions, i.e. pH values from 4 to 7.8, salt concentrations 
ranging from 0.2 M to 1.2 M and also include various additives, such as glycerol, DMSO 
and guanidine salt. This is important because tuning the solution conditions by salts and 
additives does not only affect protein-protein interactions, but also protein stability against 
aggregation and denaturatiod ^h 4 | i 5 i 99 Hi 02 ] protein adsorption il 22 Hinzl ^^d therefore, 

insights into protein-protein interactions might in turn also be relevant for these aspects. 

Published phase diagrams often give protein mass concentrations Cp [mg/mL] only. Hence, 
protein volume fractions (j) have to be calculated. For all data considered here, we have 
calculated 0 from Cp assuming 0 = Cp/pp, where Pp = 1.351 g/cw? is the protein mass 
density.l^^^S^® Note that slightly smaller values for pp are sometimes used in the literature. 
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Table I. Solution conditions under which the experimental gas-liquid binodals of lysozyme solutions 
have been determined: pH value, NaCl concentration Cg, buffer concentration Cb, additives (glycerol, 
DMSO, guandine hydrochloride). Temperatures at the critical point Tc are obtained by fits based 


on Eq. (14), typical error ATc 


from experiments. 


0.5 K. Second virial coefficient at the critical point, b^, if known 


ref. 


pH 

M 

Cb /M 

additives 

Tc/K 

^2 

24 


4.5 

0.7 

0.05 

- 

291.5 

-1.8 

24 

❖ 

4.5 

0.7 

0.05 

10% glyc. 

278.8 

-1.8 

24 

'v. 

4.5 

0.7 

0.05 

20% glyc. 

267.4 

-1.8 

24 


4.5 

0.7 

0.05 

5% DMSO 

281.4 

-1.8 

24 

- 

4.5 

0.9 

0.05 

- 

298.5 

-1.8 

24 

© 

4.5 

0.9 

0.05 

10% glyc. 

284.0 

-1.8 

24 

A 

4.5 

0.9 

0.05 

20% glyc. 

271.4 

-1.8 

52 

□ 

7.8 

0.2 

0.02 
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Figure 3. Experimental gas-liquid binodals of protein solutions (small symbols) compared to 
simulation data (large symbols), (a) Representative experimental datapl fitted by RGT scaling 
of Eq. ( [I4| ) (dashed lines) in order to estimate the critical points {(pcTc), which are indicated by 
crosses, (b) Experimental data from (a), scaled to Tc. (c) Experimental data for various solution 
conditions (Tab. [T]) and simulation data for SW fluids with A = 1.05 and 1.1 (large symbols), 
scaled to Tc. (d) Experimental binodals from (c) vs. effective volume fraction ^efj, scaled assuming 
a temperature-independent effective diameter Ueff, and SW fluid binodals that do not require 
scaling. Experimental data taken 

e.g. based on the molecular weight or different ways to quantify the volume of a lysozyme 
molecule in solution.l2^I^^*S3 

To illustrate gas-liquid binodals of protein solutions on an absolute temperature scale, 
we use the binodals of lysozyme solutions in the presence of 0.9 M NaCl as well as various 
amounts (vol.%) of glycerol (Fig. [^, symbols). In order to estimate the critical point from 
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experimental data, binodals of protein solutions are often fitted by 


r(0) = T, 




( 14 ) 


where a is a fitting parameter and f3 = 0.325 is the critical exponent for binary demixing 
from renormalization-group theory (RGT).® We note that this equation, based on critical 
phenomena, is valid only in the very proximity of the critical point, where the different shapes 
of the gas and liquid branches can be neglected. Previous experiments have suggested that 
0.15 0c ^ 0.21 (Fig. [^, colored symbols at an arbitrary short-range value of y) EIMSMIIMl 
which is slightly lower than those of the short-range SW fluids (black symbols: 0.2 ^ 0c ^ 
0.28). To ensure that only data in the very vicinity of the critical point are included, the 
htting procedure is restrained to experimental data in the region 0.05 < 0 < 0.35. For data 


at smaller or larger values of 0, the binodals show signihcant deviations from Eq. (14) which 
would systematically distort the htting. The data shown in Figure have been htted in this 
way (dashed lines) and the estimated critical points are indicated by crosses. A qualitative 
inspection reveals that the data around the critical point are well htted, but indeed, at very 
low and very high volume fractions, we observe deviations between experimental binodals 


and RGT hts. This is due to the fact that Eq. (14) assumes a certain width (0) and a 
symmetric binodal with respect to the critical point. While some studies reported symmetric 
binodals with respect to the critical pointa battening of protein binodals in the liquid 
branch (at high volume fractions) can be observed in Figure in agreement with other 

studies.™® 


The RGT hts to cloud-point measurements (Fig. |^) yield values of the critical tem¬ 
perature Tc, which can be estimated with a typical uncertainty of ATj. cs 0.5 K (Tab. |^. 
Therefore, the critical temperature turns out to be a reasonable scaling parameter. In 
contrast, it is far more difficult to estimate reliable values of the critical volume fraction 0c, 
in parts due to the occasionally observed battening of protein binodals. Typical uncertain¬ 
ties obtained from the htting procedure are A0c ~ 0.02 — 0.06, i.e., up to ~ 30 % relative 
uncertainty. Therefore, a meaningful scaling of experimental data to the critical volume 
fraction can be difficult. 

In Figure |^, we have scaled the data of Figure to the critical temperature Tc and 
indicated the critical points as well as their estimated error bars. The diherent binodals have 
a similar shape, but their gas branches do not collapse, which is in contrast to the binodals of 
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short-range SW fluids (Fig. We then have included further data sets with very different 
solution conditions (Tab.[^ and scaled their experimental binodals to their respective critical 
temperatures (Fig. [^). Again, the scaled experimental data do not collapse onto a single 
master curve; at a given volume fraction 0 > 0.05, we rather observe a certain distribution 
of T values {AT/Tc ~ 0.01). A comparison with (simulations of) SW fluids is difficult, 
because most experimental data of protein binodals are limited to a very narrow temperature 
window (0.95 < T/T^ < 1), while simulation data in this temperature window, i.e. close to 
the critical point, is scarce. Hence, to allow for a direct comparison, we added our new MC 
simulation results for SW fluids with A = 1.05 and 1.1 for T < (Fig. [^, large symbols). 
The gas-liquid binodals of SW fluids are narrower than those of protein solutions. 


Different approaches to compare SW fluids with experiments on proteins have been fol¬ 
lowed in literature due to the variety of possible scaling parameters (e.g., Tc,0c,e, A). Lo¬ 
makin et al.^ aimed at modelling liquid-liquid phase separation of protein solutions. Match¬ 
ing the volume fraction at the critical point of the protein solution, they chose a SW fluid 
with A = 1.25. However, the binodal of the corresponding SW fluid is much narrower than 
the experimentally determined one. They thus introduced a temperature dependent attrac¬ 
tion strength e(T) = ec[l -|- 7 (T — Tc)/Tc] (with 7 = —3) to “broaden” the coexistence 
curve of the SW fluid. We have transferred this approach to model the temperature depen¬ 


dence of experimentally determined second virial coefficients of protein solutions (Sec. IV 


Fig. §3, dotted line). However, this procedure works only if the values of experiments and 
simulations agree with each other (cf. Fig. [^). Grigsby et al.l^ used cloud-point temper¬ 
ature (CPT) measurements and a perturbation theory approach in order to estimate how 
the strength of the attraction e of a SW fluid with A = 1.2 varies for different salts as a 
function of their concentration. Wentzel and Guntoii^ scaled experimental protein and 
simulation binodals (A = 1.15) at a specific number density and used GPT measurements at 
this density to calculate the liquid-liquid coexistence surface in protein-salt solutions and 
to determine the salt-dependence of the second virial coefficient. However, the predicted 
coexistence surfaces seem to be too narrow compared with experimental data (as in Fig. |^) 
and this model implicitly assumes that the temperature dependence of experimental 62 val¬ 
ues is well described by the corresponding relation of SW fluids (cf. Fig. |^). As suggested 
by the van der Waals corresponding states law, often binodals of SW fluids and protein 
solutions are not only scaled to the critical temperature, but also to the critical volume 
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fraction, although the value of 0c can only be determined with a large uncertainty in 
many cases. Nevertheless, DudaP^ could show that the binodal of very short-range SW fluids 
(A = 1.05 and 1.10) nicely agrees with experimentally determined binodals of lysozyme and 
7 -crystallins (and also with a binodal of a patchy fluicP^, while binodals of SW fluids with 
A = 1.15 are too narrow and only qualitatively capture features of the protein binodals.®^ 
A qualitative difference between the SW potential and the protein-protein interaction 
potential is the presence of a repulsive contribution to the latter. According to the Noro- 
Frenkel ELCS,^it should be possible to map the binodals of protein solution onto effective 
SW fluids if the repulsive contribution of the interaction potential is taken into account in 
terms of an effective hard-core diameter aes, as defined in Eq. (|^. However, the protein- 
protein interaction potential U{r) is usually not known and therefore this procedure, i.e., 
to calculate defr from U{r), has not been applied to experimental data so far. Here, we 
determine by a comparison of experimental binodals with our simulation data of the 
SW fluid with A = 1.05. We £t each gas branch of the protein binodals (Fig. [^) to the 
simulation data using only one £t parameter, namely cTefr. This implies that aes is assumed to 
be temperature-independent, although, in general, the effective hard-core diameter might be 
temperature dependent. (In Sec. |V^ (Fig. [^), we compare the values of with calculations 
of (j^^°(T) based on DLVO theory.) The parameter cJefr might be related to the fitting 


parameter a of Eq. (14), describing the width of the binodal. However, as the value of a is 
strongly correlated with the value of 0c, which is retrieved by fitting, a is not known with 
high accuracy. 

Figure shows the scaled protein binodals as a function of effective volume fraction 
0eff = ((Jeff/o')^0. All experimental data fall on a single master curve in the gas branch, 
which - by construction - is given by the binodal of a short-range SW fluid. Since the 
binodals of the short-range SW fluids almost collapse in the low density part (Fig. [^), this 
is, to a good approximation, valid for all short-range SW fluids with A < 1.25. Interestingly, 


the assumption of a temperature-independent hard-core diameter seems to work quite well. 
We note that the experimental data in the liquid branch show a wide spread. While, for some 
data sets (blue, green and magenta symbols), we observe T/Tc > 0.995 even for 0 > 0.4, 
which corresponds to a significant flattening of the liquid branch, other data sets (red and 
dark yellow symbols) show roughly symmetric gas-liquid binodals almost agreeing with the 
liquid branch of the SW fluid. The differences of the experimental binodals at high volume 
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fractions are mainly ascribed to experimental difficnlties in the determination of CPT at 
large (j)W^ 

To test whether the experimental data obey the Noro-Frenkel predictions, they need to 
be plotted as a fnnction of second virial coefficient 62 instead of temperatnre T (Sec. 0 - 


Thns, we hrst discnss the temperatnre dependence of 62 (Sec. IV). 


IV. TEMPERATURE DEPENDENCE OF THE SECOND VIRIAL 
COEFFICIENT OF PROTEIN SOLUTIONS 

In order to present gas-Iiqnid binodals in the 62 — plane, one has to replace the tem¬ 
peratnre T with the second virial coefficient 62 , he., the relation h 2 {T) between temperatnre 
and second virial coefficient has to be known at all clond-point temperatnres. In practice, 
the second virial coefficient is measnred at a conple of temperatnres only, becanse the de¬ 
termination of 62 valnes, e.g. by static light scattering, represents an experimental challenge 
dne to the small size of proteins and the drastic effects of imparities or aggregates on the 
scattering intensity. Typical experimental nncertainties of 62 valnes are on the order of 
A 62 ~ 0.2 — 0.4.^ In Fignre]^, we show & 2 (^) measnrements corresponding to the sointion 
conditions of the binodals of Fignre as a typical example.*^ In many cases, an approxi¬ 
mate relation & 2 (T) is established based on an interpolation or even an extrapolation of the 
measnred 62 valnes. As the fnnctional form of & 2 (T) is not known a priori, henristic htting of 
lineaif^SEH or quadratic fnnctiond^to & 2 (T) data has been proposed. Moreover, it has been 
proposed that the stickiness parameter r = 1/[4(1 — 62 )] of protein solntions follows a linear 
temperatnre dependence .^2! We nse linear and qnadratic fnnctions to £t the data shown in 
Fignre illnstrating the effects of different proposed h 2 {T) relations. Depending on the 
htting procednre, systematically different 62 valnes might be predicted, in particnlar, if an 
extrapolation of the measnred data to higher and especially lower temperatnres is needed. 
The discrepancies between henristic htting approaches call for a theoretical reasoning of the 
temperatnre dependence of 62 - 

In Sec. |III[ we have shown that protein binodals lie near a master cnrve if the temperatnre 
axis is scaled to the critical temperatnre Tc (Fig. |^,d). Therefore, one might specnlate that 
the & 2 (T) curves (or h 2 {T) curves) for diherent solution conditions he near a master curve 
when the temperature axis is scaled to U,, in particular for T < U,, such that the ELCS 
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Figure 4. Temperature dependence of the second virial coefficient 62 of protein solutions, (a) 
Representative experimental data (symbols, solution conditions as in Fig. Heuristic linear and 

quadratic fits (dotted and solid lines), (b) Experimental data for various solution condition^^^EMS] 
(symbols) scaled to the critical temperature Tc of the corresponding gas-liquid binodal (Tab. 
Fig. [^), data for SW fluids with A = 1.05,1.1,1.25 (lines, top to bottom) and A = 1.25 with 
e(T) (dotted line).^ (c) Second virial coefficient of the baseline model, b^\ as obtained from the 
Donnan scaling (Eq. (d) Reduced second virial coefficient that accounts for effective 

volume fraction (/)eff (from Fig.[^). Eits to Eq. (16) assuming SW fluids with a T dependent depth 
(lines). 


mapping on equivalent SW fluids is fulfilled. In this case, experimental data scaled to the 
^2 ~ plane would fall onto the (grey-shaded) area of SW fluids in Figure 

We have collected b 2 {T) curves for a number of solution conditions in the literature ,^2312352] 
for which also the critical temperature of the gas-liquid binodal is knownP^lESl ^nd at least 
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two 62 values (per solution condition) in the vicinity of the critical point (0.95 < T/Tc < 1.05) 
are available in order to avoid extrapolation errors. Moreover, we focus on data in the high- 
salt regime (cg > 0.3 M),®^ where the range of electrostatic interactions is at most about 
20 % of the particle size and hence the interaction potential sufficiently short-range such 
that the ELCS applies. It is important to note that the experimental values of 62 at the 
critical temperature, —1.7 < &2 — are systematically lower than the values of 

short-range SW fluids, —0.8 < < —1.6 (Fig. HD- (We note that Rosenbaum et have 

also estimated values of —0.8 and —1.2 for specific solution conditions.) Nevertheless, we 
have scaled the collected h 2 {T) curves to the respective critical temperatures (Fig. |^). 

The data of Gogelein et ah®' (with the exception of one data set measured only at 
high T/Tc values) and those of Gibaud et al.^ collapse on two different master curves 
(Fig. |§D, blue and green symbols, respectively). The two data sets are based on very 
different salt concentrations and pH values. To take these differences into account, we recall 
the Donnan argument, which has been used to explain salt and charge scaling of protein 
interactionsIn particular, a scaling collapse of the crystallization boundary and the 
second virial coefficients of lysozyme solutions at room temperature has been observed 
WarrerP® has proposed a theoretical background for the charge-salt concentration scaling 
of protein interactions. He particularly assumed that the free energy of a protein solution is 
given by a baseline model, e.g., the AHS model or Sear’s sticky particle mo del,plus terms 
depending on the charge and concentration of proteins and the salt concentration. In this 
way, the second virial coefficient in the high-salt limit is given by 

= -Sf + (15) 


where R® is the second virial coefficient of the baseline model and Qe the charge of the 
protein with e the electronic charge, which is known from titration experimentsJ ^^^*^ ^ At 
room temperature (and high salt concentration), Warren found 6 ® = ^ —2.7±0.2. 

A similar scaling is found for computer simulations of the so-called primitive modeh^H^ As 
the argument for Eq. (15) also holds at another temperature, we expect another value 
and thus we assume to hnd an almost universal 6 ®(T) curve. Based on the specific solution 
conditions, namely Cs and the pH (or we have calculated bf^ values for the data 

shown in Figure]^. Indeed, in Figure]^, we observe that, at a given temperature T, data 
from different solution conditions approximately take the same b^^ value, which is - at 
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room temperature - in agreement with the value determined by Warren. (For clarity, we do 
not include further literature dataP^^®^ in this hgure, which would lead to an even larger 
spread of the data; the spread of 63 °^ is ascribed to the crude assumption of a quasi-universal 
baseline model and to experimental difficulties in determining 62 values.) The two different 
master curves observed in Figure]^ (of Gogelein et al.^ and Gibaud et al.^ are thus due 
to the very different solution conditions and follow roughly the Donnan scaling. 

The solid lines in Figure show the temperature dependence of SW fluids with A = 
1.05,1.1,1.25 (top to bottom). The temperature dependence of &2 for the short-range SW 
fluids is much weaker than that of the protein solutions. Nevertheless, one could relate 
the experimental b 2 (T) curves to SW fluids by introducing a temperature dependent range 
A(T), as suggested by Noro and Frenkel.^ It has also been proposed to introduce a linear 
temperature dependence of the SW depth e(T) (Fig. |^, dotted line),^ which was motivated 
to compensate for the too narrow width of SW binodals (compared to those of protein 
solutions). 

The interactions in SW fluids can account for attractions only, while protein interactions 
comprise attractive and repulsive contributions. To take repulsive protein-protein interac¬ 
tions into account, Noro and FrenkeP^ have proposed to introduce an effective hard-core 


diameter Ueff. We have determined values of Ueff for all protein binodals in Sec. Ill For the 
data shown in Figure]^, we compute the repulsion-corrected 62 , i-e. b^, based on Eq. ([^; 
the results are shown in Figure]^. (Note that the b^ values slightly depend on the value of 
the protein mass density pp used to calculate 0 ; a smaller pp leads to slightly smaller 63 •) 
All data sets in Figure]^, independent of their solution conditions, almost follow a master 
curve in the b^ — T/T^ representation, although the green data setd^ show slightly higher 63 
values for T < Tc and lower one for T > Tc than the blue ones.^ Since 63 takes all repulsive 
contributions into account, this includes the electrostatic effect discussed above (Eq. lUi)- 
We then propose to model the b^iT/Tc) curve using the temperature dependence of a SW 


fluid (Eq. (11)) assuming a linear T dependence of the SW depth: 


^ I = 1 - (A» - 1 ) I exp 


aT/Tc + b 

T/Tc 


- 1 


(16) 


where we set A = 1.05 and use a and b as fitting parameters. Assuming that all data points 
fall on a single curve, we fit this expression to all data (Fig. |^, solid orange line) and obtain 
a = —3.3 and b = 5.9. If it is fitted to the data sets individually, the results (dashed lines) 
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show only slightly better agreement with the data. 

It has been argued that (experimental) 62 values of protein solutions might only show a 
weak temperature dependence, mainly governed by the Boltzmann factor, and that other 
temperature effects, e.g. of hydration, are of minor importance.®*^ Despite the presumably 
weak temperature dependence of 62 , h is conceivable that, for T S> Tc, signihcant deviations 
from the master curve (Fig. |^, orange line) occur for other solution conditions. Since our 
analysis contains a limited number of data sets only, this remains to be tested in future 
experiments. 


V. THE EXTENDED LAW OF CORRESPONDING STATES FOR 
PROTEIN SOLUTIONS 


In Sec. Ill we have shown that the gas-liquid binodals of protein solutions scaled to the 
critical temperature fall on a single master curve if a temperature independent effective 
hard-core diameter aes is introduced (Fig. [^). In Sec. |IV| we have argued that the temper¬ 
ature dependence of the second virial coefficient & 2 (F) or b^iT) can be explained based on 
a SW ffuid with a T dependent depth. In view of this analysis, we are now able to present 
the experimental binodals of protein solutions in the — 4>es plane. 

Based on the master curve for the dependence of b^ on T (Fig. |^, orange line), we obtain 
a master curve of experimental protein binodals in the b^ — (j) plane (Fig. orange stars). 
If one allows for slight variations of the T dependence of b^ for different solution conditions 
(as implied by the individual hts in Fig. ii). each binodal is scaled separately (colored 
symbols), including data sets with low salt content. The data based on individual hts of 
62 (F) almost collapse onto the universal curve (orange stars), indicating that the universal 
63 (F) curve indeed provides a reasonable description of the data in the relevant temperature 
window. In addition to the experimental data in the 63 ~ 4>es representation, we re-plot the 
Noro-Frenkel area of short-range SW huids (from Fig. as well as the simulation binodals 
of the AHS huid and our SW huid data (A = 1.05,1.1) on the same axes as in Figure 
(taking into account that, in these cases, (pes = 4>)- The experimental protein data are 
completely encompassed by the (grey-shaded) Noro-Frenkel area of short-range SW huids. 
It is striking that the experimental data seems to follow very closely the binodal of the AHS 
model (except for very high volume fractions 0eff ^ 0.4). The deviations at high volume 
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Figure 5. Gas-liquid binodals, shown as vs. (j)es- comparison of experimental data (open colored 
symbols) and simulation results of the adhesive hard-sphere (AHS) fluid (black symbols) and short- 
range square-well (SW) fluids (red symbols with black frames) with the grey shaded Noro-Frenkel 
(NF) area. Data taken fronJ^ElESISSlMl^ 


fractions might be related to the appearance of non-eqnilibrium states.^ ^^ * ^^ * 


Figure [^represents the main result of this work. For a broad set of experimental condi¬ 
tions, it shows that, based on a one-parameter scaling of experimental protein data to short- 


range SW fluids, which yields aes (Sec. Ill), and an adequate description of the temperature 


dependence of the second virial coefficient, b^iT/Tc) (Sec. IV), protein phase behavior ap¬ 
pears to quantitatively agree with the predictions of the ELCS, as proposed by NF.^Thus, 
the binodals of protein solutions can be mapped onto those of short-range SW fluids if repul¬ 
sive particle-particle interactions are taken into account in terms of an effective hard-core 
diameter cJefr. Moreover, the proposed way of quantifying repulsive protein-protein inter¬ 
actions is able to clarify the discrepancies between short-range SW fluids and experimental 
protein data (Sec. [n| and Fig. m-d). 
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VI. DLVO MODEL OF PROTEIN-PROTEIN INTERACTIONS 


DLVO theory is frequently used to model protein-protein interactionsj ^^ * ^^*^^ * ^^^ ! The 
DLVO approach is more complex than SW fluid modelling, but DLVO theory contains 
attractive van der Waals as well as repulsive electrostatic interactions, which appears to be 
more realistic than simple SW fluid models. 

The DLVO interaction potential DDLVo(’")^is dehned as 


f^DLVo(^^) = Uns{r) + Usc{r) + UyBw{r). 


(17) 


The hard-sphere contribution Dhs(?’) is given by 


oo, r < a, 

Uusir) = ; 

0 , r > a. 

For r > a, the screened Coulomb contribution t/sc(’") is given by 

(Qef exp[-K(r-a)] 

- 4«„..r (l + Ka/2)= 

and the van der Waals contribution t/vDw(’") by 


(18) 


UyDw{r) = —— 


A f 


2 —2 S" 2 In 


1 - ^ 


(19) 


( 20 ) 


12 

where Cq and eg are the permittivity of vacuum and the solvent, respectively. The permittivity 
of water has a weak temperature dependence, which we describe by ■ (T — 

273.15 K)/K with Cg^o = 87.62 and Cg^i = 0.36 for 270 K < T < 330 K, in agreement 
with previous work.^^^ For water-glycerol and water-DMSO mixtures, we obtained similar 
relations.^ The Debye screening length k~^ is - for monovalent salts - given by 




e^AA / Q4>Pp 


1-0 e^eok^T \ 


-f- 2cg -1- 2ci; 


( 21 ) 


where Cg and Cb are the molar concentrations of salt and (dissociated) buffer, Na is Avo- 
gadro’s number, and the molecular weight of the protein. The strength of the van der 
Waals contribution is given by the Hamaker constant A, which for lysozyme interactions in 
brine is set to A = 8.3 k^T in agreement with previous studies.^^^^ 

It has been shown that DLVO-based models are able to quantitatively account for protein 
phase behavioi^i^l^nEll jf suitable values for the Hamaker constant are chosen. In previous 
work,^ it has been shown that DLVO-based calculations of 62 values are in very good 
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Figure 6 . Second virial coefficient 62 as a function of salt (NaCl) concentration Cs at room tem¬ 
perature: determined by scattering experiments (open symbols) DLVO calculations based on 


Eq. (22) (lines), and from cloud-point temperature (CPT) measurements (crosses, Sec. VII). 


agreement with experimental for pH 4.5 at room temperature. We now 

compare DLVO calculations with further experimental data and include data at pH 7.8. 
The second virial coefficient of the DLVO potential, can be computed directly by 


lDLVO 

"2 


= 1 + ^ 


o-" 


' ( 7-\-5 


exp 


f^sc(^) + f^VDw(^) 
knT 


r^dr. 


( 22 ) 


where 6 represents a cut-off length which is introduced to avoid divergence of the integral. 
For consistency with previous work,®E31 ^ = 0.1437 nm at room temperature T = 

293.15 K. 

In Figure]^ we compare experimental &2 values, discussed in this work (open symbols) 


with calculations (lines) based on DLVO theory (Eq. (22)). We hnd almost quantitative 


agreement between experiments and theory for both pH values and all salt concentrations 
studied. For clarity, we do not re-plot the experimental which have 

been shown to closely follow the DLVO calculations at pH 4.5.^ 

As DLVO theory is able to quantitatively account for the salt concentration and pH 
dependence of the second virial coefficient of protein solutions, we now use DLVO theory to 
determine the effective hard-core diameter in a way that is independent of our one-parameter 
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fitting approach (Sec. III). If the solution conditions are known (Tab. |^, the temperature 
dependent effective hard-core diameter within the DLVO model (j^^°(T) can be computed 
based on Eq. (§ and the repulsive part of the interaction potential Uusir) + Usc{r). We 
have calculated ct^^*^(T) in the relevant temperature window 0.9 < T/Tc < 1.2. The 
results are shown for various solution conditions in Figure]^ (lines). For comparison, we 
also show the temperature-independent values of aes, as determined by our scaling approach 
in Sec. |III| (symbols, at an arbitrary temperature). For the temperatures studied, the values 
of he in a very narrow range (1.04 < < 1.09) and show only a very weak 

temperature dependence, whereas the values of cJeff are larger than those of and lie 

in a broader range (1.16 < aefr/cx < 1.34). Nevertheless, the weak temperature dependence 
of suggests that the assumed T-independence of aes represents a reasonable hrst 


approximation. 


Having calculated the effective hard-core diameter within the DLVO model, we 

have performed a similar analysis with the experimental binodals and b 2 {T) values as in 
Sec. IV and |v| but based on ct^^°(T) instead of aes- As the values of are much 

lower than those of cxefr, the protein binodals, when scaled to the critical temperature and 
the effective volume fraction 0^^®, do not collapse with the binodals of short-range SW 
fluids, but, for a given temperature, lie at smaller 0eff (not shown). However, the (T) 

data again he near a master curve, which, following Eq. 0 . is shifted to more negative 
values (not shown). Based on these calculations, we present the experimental binodals in the 
^*DLVo _ 0 ^vo pia^j^e together with the grey shaded Noro-Frenkel ELCS area and simulation 
results of AHS and SW fluids (Fig. [^). As expected, the binodals in the 
plane are, in general, lower and toward smaller than those obtained based on cJefr 

(Fig. §. They therefore tend to he at or below the lower limit of the grey-shaded area of 
the ELCS mapping predictions. 

Thus, DLVO theory seems to consider only parts of the repulsive contributions to protein 
interactions (with respect to the ELCS mapping). A similar argument with respect to the 
repulsive interactions as quantihed by DLVO theory has been proposed by Broide et al.^^and 
others.^^^^^^ They argued that DLVO based calculations are unable to explain the stability 
of protein solutions against aggregation.® Furthermore, they suspect that repulsive protein 
hydration forces signihcantly contribute to protein-protein interactions. Whether a similar 
argument applies here needs to be examined in further experimental and theoretical studies. 
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Figure 7. (a) Effective hard-core diameter cjefr, taking into account repulsive interactions, as a 
function of temperature T: temperature-independent values determined by scaling approach (sym¬ 
bols, at an arbitrary temperature) and values determined assuming DLVO interactions (lines), 
(b) Gas-liquid binodals, shown as vs. comparison of experimental data (open 

colored symbols) with simulation results of AHS fluids (black symbols) and short-range SW flu¬ 
ids (red symbols with black frames) with the grey shaded Noro-Frenkel (NF) area. Data taken 

fron J4l24l29l52l64l 


VII. ESTIMATING SECOND VIRIAL COEFFICIENTS BASED ON 
CLOUD-POINT MEASUREMENTS 


Based on our analysis of gas-liquid binodals (Sec. Ill) and the temperature dependence 


of second virial coefficients of protein solutions (Sec. IV), we propose a method that allows 
for a quantitative estimation of 62 values at a temperature T near the critical temperature 
Tc based on cloud-point measurements. To this end, it is necessary to measure cloud point 
temperatures at various protein concentrations at the respective solution condition, so that 


the critical temperature Tc can be inferred based on Eq. (14) and the effective hard-core 


diameter Ueff can be determined by our scaling approach (Sec. III). If T is not much smaller 
or larger than Tc (0.95 < T/Tc < 1.05), possible deviations from the universal fit curve 
to b 2 {T/Tc) are expected to be small. Therefore, & 2 (T/Tc) = {aes/cr)^ b^iT/Tc) can be 
calculated if Tc and Ueff are known. Following this approach, we have determined 62 values 
at room temperature T = 293.15 K for the binodals shown in Figure]^, which only contain 
NaCl as an additive.^^^^^^ The results are shown in Figure]^ (crosses) as a function of salt 
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concentration and pH value. In the high-salt regime (cg > 0.3 M), i.e. the short-range limit, 
the results of the model calculations are in quantitative agreement with experimental data 
(symbols) as well as with the independent DLVO-based calculation of 62 values (lines). We 
have tested our approach also for the other solution conditions of the experimental 62 values 
discussed in this work (not shown) 1241521 given above, deviations between 

model calculation and experimental values are on the order of the experimental error, while 
for T S> T(. and for large values of deviations are signihcant. 

Our approach to calculate 62 values, based on CPT measurements and a comparison 
with our SW fluid data, should also apply to other solution conditions, in particular, in the 
presence of additives and in the high-salt limit, where the experimental determination of 62 
values can be difficult. 

VIII. CONCLUDING REMARKS 

In this work, we presented a direct and quantitative comparison of the experimentally 
observed binodals of protein (lysozyme) solutions and short-range square-well (SW) fluids. 
The temperature axis of the experimental binodals was scaled to the second virial coefficient 
and the concentration axis with the effective hard-core diameter which takes into account 
the repulsive contributions to the interactions. For simplicity, we assumed that the effec¬ 
tive hard-core diameter is temperature independent. It was determined by a comparison of 
experimental binodals of the protein solutions to our simulation binodals of SW fluids. Fol¬ 
lowing this procedure, the experimental data sets quantitatively agreed with the binodals of 
short-range SW fluids and of the AHS fluid, corroborating the extended law of correspond¬ 
ing states as proposed by Noro and Frenkel.lSSI Thus, our results demonstrate the predictive 
power of the ELCS also for systems as complex as protein solutions. 

If the effective hard-core diameter and hence the repulsive parts of the interactions are 
taken into account, the second virial coefficients of the protein solutions follow a master curve 
for temperatures close to the criticial temperature. Based on this hnding, we have proposed 
an approach to estimate second virial coefficients based on cloud-point measurements only. 
The obtained 62 values for various salt concentrations and pH values are in quantitative 
agreement with results from scattering experiments and DLVO-based calculations. 

Our Endings call for systematic experiments to study both the binodals of protein so- 




Table II. Monte Carlo simulation results for the gasMiquid binodal of SW fluids with A = 1.05 
(top) and A = 1.1 (bottom) for various scaled temperatures T* and densities of the gas and liquid 
phase, pq and pL, respectively. 


T* 

PG 

PL 

0.3575 

0.13393 

0.82045 

0.36 

0.16147 

0.78539 

0.3625 

0.2047 

0.7439 

0.365 

0.27549 

0.67303 

0.36679 

0.4748 

0.4748 

0.46 

0.109 

0.826 

0.465 

0.1438 

0.81346 

0.4675 

0.15931 

0.79708 

0.47 

0.17854 

0.76189 

0.4725 

0.21401 

0.73456 

0.475 

0.248 

0.687 

0.47272 

0.47961 

0.47961 


lutions and the temperature dependence of 62 for further solution conditions, in particular 
different additives. 


APPENDIX 

In Table [TT| we list our new simulation data for SW fluids with A = 1.05 and 1.1. 
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